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We present an analytic description of numerical results for the Landau-gauge SU(2) gluon prop- 
agator D{p^), obtained from lattice simulations (in the scaling region) for the largest lattice sizes 
to date, in d = 2, 3 and 4 space-time dimensions. Fits to the gluon data in 3d and in 4d show 
very good agreement with the tree-level prediction of the Refined Gribov-Zwanziger (RGZ) frame- 
work, supporting a massive behavior for D{p^) in the infrared limit. In particular, we investigate 
the propagator's pole structure and provide estimates of the dynamical mass scales that can be 
associated with dimension-two condensates in the theory. In the 2d case, fitting the data requires 
a non-integer power of the momentum p in the numerator of the expression for D{p^). In this 
case, an infinite-volume-limit extrapolation gives -D(O) = 0. Our analysis suggests that this result is 
related to a particular symmetry in the complex-pole structure of the propagator and not to purely 
imaginary poles, as would be expected in the original Gribov-Zwanziger scenario. 
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I. INTRODUCTION 

High-precision data from lattice simulations are a key ingredient in our understanding of the low-energy aspects of 
Yang-Mills theories associated with color confinement. In fact, whereas new insight into the confinement mechanism 
may be gained by investigating the properties of gauge-field configurations produced in the simulations (see e.g. [1]), 
specific features of proposed confinement scenarios may be tested by comparison with lattice data. In this case, one 
may obtain physical values for a model's parameters by fitting the predicted expression of a given observable to its 
numerical realization. One may also hope to over-constrain the proposed analytic forms, if the fits can be done with 
sufiiciently high number and wide range of data points, from which systematic errors have been consistently eliminated. 
In particular, this applies to predictions for the infrared behavior of gluon and ghost propagators, formulated in Landau 
gauge for S\J{Nc) gauge theory. Here we perform a series of fits to the gluon propagator D{p-^) and tost the predictions 
of the so-called Refined Gribov-Zwanziger (RGZ) framework, which differs from the scenario originally proposed by 
Gribov [2] and Zwanziger [3] through the introduction of dimension-two condensates, associated with dynamical mass 
generation [4]. Our analysis is done for pure SU(2) gauge theory. The data have been produced previously and 
discussed in [5-7] (see also [8]), but they have not been systematically fitted until now. A companion paper with 
similar fits for the ghost propagator is imder way [9] . We note that an alternative comparison of these data to analytic 
predictions was recently presented in [10] . 

The Gribov-Zwanziger confinement scenario is based on restricting the functional integration to the first Gribov 
region SI delimited by the first Gribov horizon dfl, where the smallest nonzero (and positive) eigenvalue of the Faddeev- 
Popov matrix A4 goes to zero [2, 3]. Let us recall that limiting the gauge configurations to this region was an attempt 
— made by Gribov in Rcf. [2] — to fix the gauge completely, getting rid of spurious gauge copies, known thereafter 
as Gribov copies. Now, f2 is a convex region of very high dimensionality and therefore, as the infinite-volume limit 
is approached, the increase in entropy should favor [11] gauge configurations on the surface dfl. This in turn can 
cause the infrared enhancement of the ghost propagator (which is related to the inverse of A^), inducing long-range 
effects in the theory. Indeed, in Coulomb gauge, the restriction to the first Gribov region causes the appearance 
of a confining color-Coulomb potential [2, 12]. Thus, in this scenario, formulated for momentum-space propagators, 
the long-range features needed to explain the color-confinement mechanism are manifest in the ghost propagator, 
whereas the momentum-space gluon propagator D{p^) should be suppressed in the infrared limit. Such a suppression 
is associated with violation of spectral positivity, which is commonly interpreted as gluon confinement [3, 13, 14]. In 
particular, -D(O) is originally expected to be zero [2, 3], corresponding to maximal violation of spectral positivity. The 
parametrization of this behavior as a propagator having a pair of poles with purely imaginary masses has arisen in 
the Gribov-Zwanziger approach [2, 3], in connection with the study of gauge copies. 

Lattice studies (see [8] for a recent review) have confirmed the suppression of the gluon propagator in the infrared 
limit and the enhancement of the ghost propagator at intermediate momenta. However, considering lattice sizes large 
enough to allow the investigation of the deep infrared regime, it is clear that the results of the simulations are not 
compatible with the scenario described above. Indeed — in space-time dimension d = 3, 4 — the gluon propagator 
shows a finite value as the momentum is taken to zero and the enhancement of the ghost propagator is lost in this 
limit. We note the very large lattice sizes employed in order to observe such a behavior, _L w 20 fm and larger 
[5-7, 15-18]. In any case, violation of refiection positivity for the real-space gluon propagator (see e.g. [19]) is clearly 
observed in the data. 

Recently, the quantitative description of the massive behavior for the gluon propagator has been studied by several 
groups [10, 18, 20-30], based on different proposed analytic forms. Earlier attempts of fitting gluon- propagator data 
can be found for example in Refs. [31-34]. We note that some of these studies (see e.g. [31, 34]) have considered 
the so-called Gribov-Stingl form [35, 36] for modeling the massive behavior of the gluon propagator. This form is a 
generalization of the Gribov propagator described above, including pairs of complex-conjugate poles with a nonzero 
real part. As illustrated below, the behavior predicted for D(p^) in the RGZ framework is also based on general 
complex-conjugate poles for the (massive) propagator. This proposed form is given for SU(A'^c) gauge theory and 
for four or three space-time dimensions [4, 37-40]. On the contrary, in the 2d case, the RGZ approach cannot be 
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implemented, since the dimension-two condensates would induce severe infrared singularities, precluding the restriction 
of the functional integration to the first Gribov region [41]. By fitting rational functions of to the whole range 
of data for the SU(2) gluon propagator, we are able to obtain estimates for the physical values of the masses in the 
RGZ framework, as well as to gain a better understanding of the pole structure in the proposed expressions. In each 
case, wc look for the best fit to the data, with the smallest number of independent parameters, and relate them to 
the condensates in the proposed analytic forms only at the end. Put differently, the predicted dependence of the fit 
parameters on the condensates is not imposed in the fitting form, but is obtained as a result of the fit. This allows 
us to use a wide fitting range, considering all data points. We note that predictions from the RGZ framework were 
already tested in [23] , showing good fits (using a somewhat different analytic form and a smaller fitting range) to 4d 
lattice data for the SU(3) case. 

The paper is organized as follows. The Gribov- Zwanziger scenario is briefly reviewed in Section II. The introduction 
of condensates as part of the RGZ scenario is summarized is Section III, where we present the expressions to be fitted 
to the lattice data. The numerical results are discussed in general in Section IV and in particular for the 4d, 3d and 
2d cases respectively in Sections V, VI and VII. We present our conclusions in Section VIII. 



II. THE GRIBOV-ZWANZIGER ACTION 



The Gribov- Zwanziger (GZ) action, introduced in 1989 [13], implements an all-order restriction of the path integral 
to the first Gribov region 



Q= {Al{x): d^Al{x)=0,M''\x,y)>Q] , 
where A'^{x) is the gauge field and j/) is the Landau-gauge Faddcev-Popov operator 

M''\x,y) = -6{x-y)d^Df = 5{x - y) {-5'^'' dl + /""^S^^^) . 



(1) 



(2) 



By introducing auxiliary fields a pair of complex-conjugate bosonic fields (^JJ'^, '^f!^) and a pair of anticommuting 
complex-conjugate fields oo^f^ — one is able to obtain a local renormalizable action [11, 42, 43]. More precisely, 

the generating functional for the GZ action can be written in d space-time dimensions as [13, 43, 44] 



Z{J) = y[d$]e^°^+/'i'^-^M(^)^J 
where ^gz is the local GZ action given by 



(x) 



Sg2 



So 



with 



and 



So = Sym + Ses + 



T^oco Ttah^.hc r-roc^ I Tif^h , ,hc\ „ I rT<^^\ £abc Tibm , ^cn 



-7 



/" ^abc Aa,J>c i £ahc Aa-r^hc i ^ / \r2 i\ ^2 

g J d X f A^ip^ + f A^ipi^ + - (W^ - Ij 7 



(3) 



(4) 



(5) 



(6) 



Here, a, b, c, m and n are color indices in the adjoint representation, Nc is the number of colors, 7 is the so-called 
Gribov parameter, Sym is the classical Yang-Mills action 



Sym ^ - I d^xF,? F,?. 



and 5gf is the Landau-gauge-fixing action 

5gf = j d^x {b'^d^A; + ^d^Dl 



(7) 



(8) 
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where the auxihary field 6° is a Lagrange multipHer enforcing Landau gauge and (c", c") are the Faddeev-Popov ghost 
fields. Also, we indicate with [d$] the integration over aU fields $ S {A^, c"*, c", 6", (y?^'^, o;"^, uJj^'^}. Notice that 
one can simplify the notation of the auxiliary fields {^p*^, tp^f, uJ^f, w^") in the action Sq using the symmetry of this 
action with respect to the composite index i = {fi, c). Thus, we can set 

t-r^ac ,^ac T^ac , ,ac\ /Trra 7-ra , ,a\ /^^ 

VP IX ^'Pp. ■•'^p. ■>^ix) = iVi, fi, i^i, i^i) (9) 

and write 

So = Sym + % + J d^x [iftd^ (Dt'^'i) - oJtd. (Dfojl) - g {d^oJ^) r"^ D'J^c"^^^] . (10) 
Finally, the parameter 7 in Eq. (6) is fixed by the so-called gap equation (also known as the horizon condition) 

( gr'^A;^'; ) + ( gr'^A^ip"; ) + 2j^d {N^ - 1) = , (11) 

where ( ) indicates the expectation value in the measure defined by Eq. (3). This condition is a consequence of the 
restriction of the path integral to the first Gribov region. 

As mentioned in the Introduction, one of the main outcomes of the GZ theory is the modification of the behavior 
of gluon and ghost propagators in the infrared (IR) limit in comparison with the perturbative behavior [2, 3, 11, 
13, 14, 45]. Indeed, the gluon propagator becomes IR-suppressed with a tree-level behavior given by 

( A;{p)Ati-p) ) . Dip-) ( - ^ ) = ^^^rj^, ( - ^ ) • (12) 

This result is confirmed by one- loop calculations [46, 47]. The above expression for the gluon propagator implies that 
D{p^) is null at zero momentum, which in turn indicates maximal violation of refiection positivity for the real-space 
gluon propagator D{x) [45]. This violation is usually considered a manifestation of gluon confinement [3, 13, 14]. At 
the same time, one finds that the ghost propagator displays an enhanced IR behavior 

{c'^{p)^{p)) = S'^'Qip-) S'^"^. (13) 

This behavior is indicative of a long-range interaction in the theory and it should be related to quark confinement 
[2, 14, 48]. 



III. THE REFINED GRIBOV-ZWANZIGER FRAMEWORK 

More recently, the GZ action has been "refined" by taking into account the possible existence of dimension-two 
condensates [4, 37-40]. In the most general case [40], four different condensates are considered, i.e. 

{A;A;)^-m' {ipl^D^M- {^pI^D^p (14) 

where we have listed the dynamical mass associated to each condensate. (Note that p is complex, whereas — and 
are real and positive.) The condensate — is directly related to the gluon condensate {g'^A^) (see e.g. [23]). 
One can show [40] that the Refined Gribov- Zwanziger (RGZ) action can be renormalized. At the same time, there 
is clear evidence that the original GZ theory dynamically transforms into the refined theory, since the minimum of 
the associated effective potential favors non- vanishing condensates [40]. As displayed below, a non-zero value for 
these condensates has an effect on the IR behavior of gluon propagators. The effect on the ghost propagator will be 
discussed in a forthcoming work [9]. 



A. The Gluon Propagator 

In the presence of the four condensates considered in Eq. (14), the GZ gluon propagator (12) is modified [40] as 

p6 + p4 (^^2 + 2M2) + p2 (2m2M2 + Af4 + A4 - ppt) + (M4 - ppt) + M2A4 - ^ (p + pt) ' 

where the condensates m^, Af^, p arc described above and A'* is related to the Gribov parameter 7 through = 
2g^Nc'f'^. Since p and are complex-conjugate quantities, we can set 

p = pi+ ip2 

= pi- ip2 (16) 

and rewrite Eq. (15) as 

2. ^ p4 + 2MV + M4-(p2 + p2) 

^ p6 + pi (^2 + 2M2) + p2 [2m2M2 + M4 + A4 - (p2 + p2)] ^ ^2 [J\^4 _ (^2 ^ ^2^] ^ ^4 (^2 _ ' ^ 

It is interesting to notice that this propagator gets simplified if p = p^ = pi (i.e. p2 = 0), which corresponds to the 
equality {'<fi(p) = ( (pip ) from (14). Indeed, in this case one can factorize the quantity + — pi in the numerator 
and in the denominator of the above formula, obtaining 

' + p2 (M2 + m2 + Pi) + m2 (M2 + p^) + A* " ^ ' 

Clearly, both propagators (17) and (18) have, in principle, a finite nonzero value at zero momentum. Nevertheless, if 
the value of £'(0) is sufficiently small, one still finds that the real-space propagator D{x) becomes negative for some 
(large) value of x, i.e. reflection positivity can also be violated for these propagators. 

Note that both Eqs. (17) and (18) can be decomposed as sums of propagators of the type a/ (p'^+co^). In particular, 
we can write Eq. (17) as 

£)(p2) = + , ^ „ + , ^ „ . (19) 

^ ' p^ + Lof p^ +LU^ p^ ^ ^ 

To this end, we only need to solve the cubic cqiiation 

a;3 + (m^ + 2M^) + x [2m^M^ + + - {pi + pi)] + [M" - (p? + pj)] + A^ (M^ - pi) = , (20) 

obtained by setting p^ = x in the denominator of Eq. (17), and to find its three roots wf, ^"^^ ^3- -^t the same 
time, the gluon propagator in Eq. (18) can be written as 

Dip'') = + , (21) 

where wc expect to have a_ = ol\ if uPl. = (w+)*, i.e. if uj\ and iJ^_ are complex conjugates. Here, uj\ are the roots 
of the quadratic equation 

a;2 + a; (M^ + + pi) + (M^ + p^) + A* = , (22) 

obtained by setting j? = x'vs^ the denominator of Eq. (18). Clearly, one finds complex-conjugate poles if \M'^ — rr? -\- 
P\\ < 2A2. 

Let us remark that rational forms such as (17) and (18) for the gluon propagator were considered by Stingl [35, 36], 
as a way of accounting for nonperturbative effects in an extended perturbative approach to Euclidean QCD. More 

precisely, in his treatment, one expresses the proper vertices of the theory as an iterative sequence of functions yielding 
a self-consistent solution to the Dyson- Schwinger equations. In particular, for the gluon propagator, this sequence is 
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written [sec Eq. (2.10) in Ref. [36]] in terms of ratios of polynomials in the variable p^, of degree r in the numerator 
and r + 1 in the denominator, with r = 0,1,2,... . This functional form is then related, via operator product 
expansion, to the possible existence of vacuum condensates of dimension 2n, with n > 1. At the same time, the 
associated complex-conjugate poles^ are interpreted as short-lived elementary excitations of the gluon field [3, 35, 36]. 

By comparison, in the RGZ framework, one proposes specific forms for the dimension-two condensates — related to 
the auxiliary fields of the GZ action and then obtains (at tree level) the rational functions in Eqs. (17) and (18), 
which correspond respectively to cases with r = 3 and 2 in Stingl's iterative sequence. 

In Section V below, we show that the simplest rational form [with r = 2, corresponding to Eq. (18)] works well in 
the 4d SU(2) case. A similar result was obtained for the SU(3) case in [23].^ It may be noted that, in Ref. [31], lattice 
data for the 4d SU(3) Landau-gauge gluon propagator were fitted using the above sequence of functions with r = 2, 4 
and it was found that a good description of the data can be achieved only for r = 4. Note, however, that the fit was 
performed for the real-space propagator, for which the analysis is known to be complicated by several technical issues 
(see e.g. [19, 51]). Moreover, although the lattice volume considered was rather large, the study employed asymmetric 
lattices, which may give rise to systematic effects [52] . 

IV. NUMERICAL SIMULATIONS 

The data presented here for the SU(2) Landau-gauge gluon propagator were produced in 2007. The 3d and 4d 
cases were run on the 4.5 Tflops IBM supercomputer at LCCA USP [53], whereas the 2d case was run on various 
PC clusters at the IFSC-USP. Most of these data have already been discussed in Refs. [5-8], but they were not 
systematically fitted up to now. 

In the 4d case, we have considered lattice sides N = 48, 56, 64, 80, 96 and 128, with lattice parameter (3 = 2.2. 
The corresponding lattice spacing a is approximately 0.210 fermi, implying that the smallest non-zero momentum 
Pmin = 2sm{TT/N) is about 46 MeV in physical units for the N = 128 lattice. In this case, the physical lattice volume 
V = N'^ is about (27 fermi)^. The number of gauge-field configurations produced was 168 for N = 128 and about 250 
for the other lattice sizes. 

In 3d, we have N = 140, 200, 240 and 320 at /3 = 3.0, with a w 0.268 fermi. Then, for the lattice volume 320^ the 
smallest non-zero momentum Pmin is about 14 MeV and the physical volume corresponds to about (85 fermi)^. The 
number of configurations was 630, 525, 350 and 125, respectively for the four lattice sizes. 

Finally, in the two-dimensional case, we considered N = 80, 120, 160, 200, 240, 280 and 320 at /3 = 10.0. In this 
case the lattice spacing is about 0.219 fermi, the lattice volumes 320^ correspond to V (70 fermi)^ and in this case 
Pmin ~ 18 MeV. We have about 600 configurations for each lattice volume. 

In all cases we set the lattice spacing a by considering the input value a^/^ = 0.44 GeV for the string tension, which 
is a typical value for this quantity in the 4d SU(3) case. The evaluation of the lattice string tension is described in 
[54], [34] and [55], respectively for d = 2, 3 and 4. Note that all our runs are in the scaling region [34, 55, 56] and all 
data refer to the SU(2) case. Possible systematic effects due to Gribov copies [57-61] as well as unquenching effects 
[62-65] were not considered here. Finite-volume effects, on the other hand, are well under control. In particular, in 
3d and in 4d, our largest lattice volumes can be already considered as infinite. In the 2d case, a simple extrapolation 
to infinite volume needs to be considered, as done in Section VII to obtain the limiting value of -D(O). 

Configurations have been generated by alternating heat-bath updates of the link variables with micro-canonical 
steps, in order to reduce the problem of critical slowing-down (see for example [66] and references therein). Gauge- 
fixing to Landau gauge was done using the stochastic overrelaxation algorithm [54, 66]. Let us also recall [see Eq. 



^ See [36, 49, 50] for some considerations concerning the issue of causality for propagators with complex poles. 

^ Let us mention that the condensate p was not considered in Ref. [23]. Therefore, when discussing fit results using Eq. (18), we must 
compare their values for to our values for + pi . 



7 



(12)] that the gluon propagator D{p'^) in Landau gauge is evaluated using 

f,-2mp-{x-v)/N / T) T) \ 

= 12- V = ^'^ [a^-^ - ^j^ip")- (23) 

x,y \ P / 

Here is the lattice gluon field defined as"^ 

where Ufi{x) are the usual link variables of the Wilson action. Also, the momentum components P/j, are given by 

p, = 2 sin (25) 

and p^ takes values 0, 1, — 1. 

In 2d we considered momenta with components (j>, 0) and (j>,p), plus all possible permutations of the components. 
Similarly, in 3d, we have data for momenta with components {p, 0, 0), {p,p, 0) and {p,p,p) and all possible permutations 
of components. Finally, in 4d, we evaluated the propagator for momenta with components (p, 0,0,0), {p,p, 0,0), 
{p,p,p,0) and (p,p,p,p). In this case, we considered all possible permutations of the components for momenta 
of the type (p, 0,0, 0). On the contrary, we did not consider permutations for the momenta {p,p,p,0) and in the 
case {p,p, 0,0) we allowed all permutations satisfying the constraint p4 = 0. When permutations of the momentum 
components were available, an average over the different permutations was taken for each configuration. In order to 
reduce discretization effects — and in particular those related to the breaking of rotational symmetry [32, 68, 69] — 
we have considered, in addition to the usual (unimproved) momentum defined by the squared magnitude of the lattice 
momenta 

p' = EpI' (26) 
ft 

the improved definition [68] 

p' = Epl + r2l2pt- (27) 

ft n 

This definition does not affect the value of p^ in the IR limit, but modifies its value significantly for large momenta. 
In particular, the largest value of p^ — obtained when = A^/2 in Eq. (25) for all directions /i — is given (in lattice 
units, for the d-dimensional case), by 4d if the unimproved definition is considered, and by 16d/3 in the improved 
case. For the j3 values considered here, this implies that the largest momentum Pmax is about 2.54, 2.55 and 3.75 GeV, 
respectively in 2d, 3d and 4d in the unimproved case, and about 2.94, 2.94 and 4.33 GeV using improved momenta. 

In the next sections we present fits (obtained using gnuplot) of the 4d, 3d and 2d data for the SU(2) gluon propagator 
and compare the fit results to the predictions of the RGZ action, discussed above in Section III A. We remark that the 
shown data for D{p'^) are not normalized. Note that a renormalization condition at a given scale /^^ would correspond 
to a rescaling of the overall factor C in the fitting forms considered below. The condensates and the poles, on the other 
hand, are not affected by such a renormalization. We also note that, whenever possible, we avoid rounding off the 
values of the fit parameters. On the contrary, values for the associated physical quantities (i.e. condensates and poles) 
are rounded to show errors with one significant digit only. We refer to the gnuplot documentation [70, "Statistical 
Overview" section] for information on the significance of the standard errors calculated for the fit parameters. 



^ With this definition of the lattice gluon field, the gluon propagator evaluated on the lattice corresponds to the propagator g^D{p'^) in 
the continuum, which has mass dimension 2— d in the d-dimensional case [67]. 
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FIG. 1. Plot of the 4d gluon propagator D{p^) (in GeV^^) as a function of the (unimproved) momentum p (in GeV) for the 
lattice volume V = 128*. As a consequence of the breaking of rotational invariance, the data do not produce a smooth curve. 
Note the logarithmic scale on the y axis. 



V. THE 4D CASE 

As a first attempt in the 4d case, we consider a fitting function of the simplest Gribov-Stingl form"^ 

Mp') = C , (28) 

which corresponds to the RGZ propagator in Eq. (18), modulo a global rescaling factor C. We note that, in order 
to improve the stability of the fit, we impose some parameters to be positive, by setting them to be squares. The 
results of the fit for all lattice volumes, using unimproved and improved momenta, are reported respectively in Tables 
I and II. From the x^/d.o.f. values one clearly concludes that the use of improved momenta makes the behavior of the 
gluon propagator smoother, allowing a better fit to the data. This is also seen by comparing the data in Figs. 1 and 
2, plotted respectively for unimproved and improved momenta. Let us stress that we are fitting the whole momentum 
range available and that, for the largest lattice volume, we have 257 data points. 

In order to extract the value of the condensates described in Section III A above, we now consider only the fit 
results for the volume V = 128^ (using improved momenta), reported in the last row of Table II. The corresponding 
plot is shown in Fig. 2. By setting fi{p^) equal to the RGZ propagator in Eq. (18) (modulo the global factor C) and 
using propagation of error, we find for the condensates the values reported in the first column of Table III. (Note 
that, for this fitting form, the condensates A'P and pi cannot be determined separately.) Let us mention that the 
values obtained here for + pi, and are in good quantitative agreement with the corresponding values — 



* Since our largest momentum is of the order of 4 GeV, ultraviolet logarithmic corrections are not important to describe the lattice data 
and they are not included in the fitting functions proposed here. This also avoids the problem of having to regularize the corresponding 
Landau pole by hand. 
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V 


C 


«(GeV) 


t(GeV2) 


s(GeV2) 


xVd.o.f. 


48* 


0.567 (0.016) 


0.507 (0.065) 


0.607 (0.023) 


2.417 (0.269) 


16.47 


56* 


0.572 (0.014) 


0.495 (0.058) 


0.602 (0.021) 


2.344 (0.234) 


15.46 


o4 


r\ rrcc f f\ ri"l 0^ 

O.OOD (U.Old) 


c\ COO f c\ rtco\ 

0.522 (0.052) 


r\ z?-io f c\ ni r\\ 

0.612 (0.019) 


2.452 (U.223) 


15.00 


80* 


0.562 (0.012) 


0.496 (0.048) 


0.612 (0.017) 


2.443 (0.199) 


16.59 


96* 


0.567 (0.011) 


0.484 (0.044) 


0.604 (0.016) 


2.367 (0.178) 


15.96 


128* 


0.560 (0.010) 


0.534 (0.037) 


0.621 (0.014) 


2.553 (0.166) 


10.65 



TABLE I. Fits of the gluon-propagator data in the 4d case, for different lattice volumes, using the fitting function fi{p^) in 
Eq. (28) and unimproved momenta [see Eq. (26)]. We report, besides the value of the fit parameters, the x^/d-o.f. obtained in 
each case. The whole range of momenta was considered for the fit. Errors shown in parentheses correspond to one standard 
deviation. 



V 


C 


m(GcV) 


f(GcV') 


s(GcV") 


xVd.o.f. 


48* 


0.791 (0.007) 


0.755 (0.027) 


0.707 (0.013) 


2.419 (0.119) 


2.09 


56* 


0.801 (0.006) 


0.734 (0.023) 


0.696 (0.012) 


2.305 (0.100) 


1.92 


64* 


0.791 (0.007) 


0.760 (0.024) 


0.710 (0.012) 


2.425 (0.108) 


2.35 


80* 


0.785 (0.005) 


0.734 (0.019) 


0.708 (0.009) 


2.404 (0.084) 


2.04 


96* 


0.795 (0.004) 


0.717 (0.016) 


0.694 (0.008) 


2.291 (0.068) 


1.66 


128* 


0.784 (0.005) 


0.768 (0.017) 


0.720 (0.009) 


2.508 (0.078) 


1.63 



TABLE II. Fits of the gluon-propagator data in the 4d case, for different lattice volumes, using the fitting function fi{p^) in 
Eq. (28) and improved momenta [see Eq. (27)]. We report, besides the value of the fit parameters, the x^/d.o.f. obtained in 
each case. The whole range of momenta was considered for the fit. Errors shown in parentheses correspond to one standard 
deviation. 

respectively indicated with M^, and 2g'^N'y'^ — reported in Ref. [23] for the SU(3) case.^ Also, as remarked above, 
the condensate may be used to obtain a value for the gluon condensate {g^A^), through the relation (see e.g. [23]) 

In our case, the value rri^ = —1.92(9) from Table III (using propagation of error) yields {g'^A'^) = 1.99(9) GeV^. 

Furthermore, we verify from Table III that jM^ — mP + pij < 2X^, justifying our expectation (see end of Section 
III A) that the propagator may be decomposed in terms of a pair of complex- conjugate poles. We can thus write [see 
Eq. (21)] 

2x "+ , a- 2ap^ + 2{av + bw) 

with a± = a±ib and = vzL iw. The results for the parameters a, b, v and w (again using propagation of error) are 
shown in the first column of Table IV. Thus, the poles are complex conjugates whose imaginary part is more than 
twice their real part. We recall that a Gribov propagator would have a null real part. 

In order to have better control of the errors on the values of the condensates and poles, we have re-done the estimates 
described above using a Monte Carlo error analysis (with 10000 samples).^ The corresponding results are reported in 
the second columns of Table III and Table IV. Finally, we repeated the fit and the evaluation of the condensates and 



° For comparison with our values in Tabic III, the SU(3) condensates from [23] are respectively 2.15(13) GcV^, -1.81(14) GcV^ and 
4.16(38) GeV*. 

^ To this end, we considered independent Gaussian distributions for the fit parameters. Thus, this Monte Carlo analysis may be considered 
as a numerical check of the analytic propagation of error. 
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FIG. 2. Plot of the 4d gluon propagator D{p^) (in GeV~^) as a function of the (improved) momentum p (in GeV) for the 
lattice volume V = 128*. We also show the fitting function /i(p^) [see Eq. (28)] with the parameters reported in the last row 
of Table II. Note the logarithmic scale on the y axis. 



parameter 


propagation of error 


Monte Carlo analysis 


bootstrap analysis 


+ pi (GeV^) 


2.51(8) 


2.51(8) 


2.3(3) 


(GeV^) 


-1.92(9) 


-1.92(9) 


-1.7(2) 


A'* (GeV^) 


5.3(9) 


5.3(4) 


4.5(9) 



TABLE III. Estimates of the parameters of the RGZ gluon propagator in Eq. (18) from fits (see last row of Table II above) to 
the equivalent form fi{p^) in Eq. (28), using propagation of error. For comparison, we also report a Monte Carlo error analysis 
with 10000 samples and a bootstrap analysis (fit results given in the text) with 500 samples. In all cases we considered the 
volume V = 128* and improved momenta. Errors shown in parentheses correspond to one standard deviation. 



poles using a bootstrap analysis (with 500 samples). In this case, from the fit of fi{p^) [see Eq. (28)] at V = 128^, 
we find the parameters 

C = 0.762 ±0.024 (31) 

u = 0.755 ± 0.035 GeV (32) 

t = 0.698 ± 0.027 GeV^ (33) 

s = 2.292 ± 0.253 GeV^ . (34) 

The corresponding results for the condensates and poles of the propagator are shown in the third columns of Tables 
III and IV. Clearly, all results obtained agree within one standard deviation. 

As a second test, we have also tried to allow for the more general form of the propagator, given in Eq. (17). To 
this end, we consider the fitting function 

+ cp^ + dp-^ + 
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parameter 


propagation of error 


Monte Carlo analysis 


bootstrap analysis 


a 


u.692[6) 


0.392(2) 


0.38(1) 


b 


1.32(7) 


1.32(5) 


1.20(7) 


v{GeY^) 


0.29(2) 


0.29(2) 


0.29(3) 


w (GeV^) 


0.66(2) 


0.66(1) 


0.64(2) 



TABLE IV. Estimates of the parameters of the function /2(p^) [see Eq. (30)] from fits (sec last row of Table II above) to the 
equivalent form /i(p^) in Eq. (28), using propagation of error. For comparison, we also report a Monte Carlo error analysis 
with 10000 samples and a bootstrap analysis (fit results given in the text) with 500 samples. In all cases we considered the 
volume V = 128* and improved momenta. Errors shown in parentheses correspond to one standard deviation. 



V 


C 


d(GeV*) 


e(GoV^) 


6(GcV*) 


c(GcV2) 


a(GeV) 


xVd.o.f. 


48* 


0.889 (0.087) 


7.742 (9.663) 


2.469 (1.651) 


26.613 (32.940) 


9.692 (13.240) 


2.085 (1.143) 


1.95 


56* 


0.798 (0.007) 


0.495 (0.020) 


0.055 (0.045) 


0.014 (0.023) 


0.581 (0.057) 


1.093 (0.030) 


2.04 


64* 


0.795 (0.010) 


0.625 (0.662) 


0.336 (0.868) 


0.543 (2.790) 


0.733 (1.031) 


1.121 (0.209) 


2.52 


80* 


0.781 (0.006) 


0.514 (0.016) 


0.059 (0.037) 


0.016 (0.020) 


0.593 (0.046) 


1.122 (0.025) 


2.08 


96* 


0.893 (0.104) 


10.539 (14.750) 


3.036 (2.233) 


39.289 (53.150) 


16.558 (24.490) 


2.720 (1.693) 


1.56 


128* 


0.784 (0.006) 


0.578 (2.192) 


0.229 (4.228) 


0.253 (9.335) 


0.691 (3.795) 


1.143 (0.843) 


1.69 



TABLE V. Fits of the gluon-propagator data in the 4d case, for different lattice volumes, using the fitting function fsip^) in 

Eq. (35) and improved momenta [see Eq. (27)]. We report, besides the value of the fit parameters, the x^/d.o.f. obtained in 
each case. The whole range of momenta was considered for the fit. Errors shown in parentheses correspond to one standard 
deviation. 



which has six parameters. As can be seen in Table V, the values of x^/d.o.f. do not improve in comparison with the 
previous (4-parameter) fit and, with the exception of the global factor C, most of the parameters are determined with 
very large errors. This suggests that the above function has too many (redundant) parameters, making the fitting 
procedure quite unstable. Next, we reduce the number of parameters by one and introduce a general form that will 
prove useful in the description of the 3d data, in Section VI. More precisely, we test the fitting function 



V 


C 


w(GeV) 


t{GcY^) 


s{GeV^) 


k{GeV^) 


xVd.o.f. 


48* 


0.802 (0.009) 


0.686 (0.081) 


0.792 (0.030) 


1.662 (0.368) 


0.547 (0.149) 


2.02 


56* 


0.809 (0.008) 


0.694 (0.063) 


0.761 (0.034) 


1.714 (0.365) 


0.622 (0.177) 


1.89 


64* 


0.802 (0.008) 


0.701 (0.071) 


0.790 (0.028) 


1.716 (0.346) 


0.573 (0.144) 


2.28 


80* 


0.793 (0.007) 


0.703 (0.046) 


0.767 (0.031) 


1.882 (0.308) 


0.668 (0.154) 


2.01 


96* 


0.804 (0.006) 


0.673 (0.043) 


0.757 (0.022) 


1.694 (0.246) 


0.625 (0.119) 


1.62 


128* 


0.793 (0.006) 


0.727 (0.045) 


0.791 (0.025) 


1.903 (0.265) 


0.631 (0.120) 


1.60 



TABLE VI. Fits of the gluon-propagator data in the 4d case, for different lattice volumes, using the fitting function fi{p^) in 
Eq. (36) and improved momenta [see Eq. (27)]. We report, besides the value of the fit parameters, the x^/d.o.f. obtained in 
each case. The whole range of momenta was considered for the fit. Errors shown in parentheses correspond to one standard 
deviation. 



J^yf > ^ (p4 + ^2 p2 + ^2) (p2 + fc) ^ p6 + (jt + U2)p4 + (fcy2 + ^2)^2 + ' ^"^"^ 

This function is of the type (35) (with different parameters), but is written as a simple generalization of fi{p^) in 
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V 


C(GeV) 


«(GeV) 


^(GeV2) 


s(GeV2) 


xVd.o.f. 


140"' 


r\ A A -1 ( r\ C\C\o\ 

U.441 (U.UUJJ 


u.SUb (u.UlS) 


v. Sao (U.UU9) 


U.217 (U.013) 


10.86 


200^ 


0.440 (0.002) 


0.305 (0.011) 


0.389 (0.008) 


0.223 (0.011) 


8.68 


240^ 


0.443 (0.002) 


0.307 (0.010) 


0.374 (0.007) 


0.198 (0.009) 


6.53 


320^ 


0.445 (0.002) 


0.296 (0.011) 


0.365 (0.006) 


0.183 (0.008) 


3.19 



TABLE VII. Fits of the gluon-propagator data in the 3d case, for different lattice volumes, using the fitting function /i(p^) in 
Eq. (28) and improved momenta [see Eq. (27)]. We report, besides the value of the fit parameters, the x^/d.o.f. obtained in 
each case. The whole range of momenta was considered for the fit. Errors shown in parentheses correspond to one standard 
deviation. 



Eq. (28).^ Fit results arc shoviTi in Tabic VI. In this case the fits look reasonable. The corresponding values for the 
condensates and A"* in Eq. (17) are obtained by a Monte Carlo analysis (with 10000 samples) using the data in the 
last row of Table VI. We note that, in this case, the fitting form allows us to evaluate M^, pi and \p\ separately. We 
find the values 

= 1.5 ± 0.1 GeV^ (37) 
= -1.7±0.3GeV^ (38) 
A* = 4.1 ± 1.0 GeV^ (39) 
pi = 0.5 ± 0.1 GeV^ (40) 
pI+pI = 0.2 ± 0.1 GeV^ . (41) 

We see that the errors are larger, and that the value of + pi is incompatible with the numbers in Table III 

(obtained assmriing p2 = 0). Also, a comparison of the values in (40) and (41) suggests a very small (and imaginary) 
value for p2, implying that p is real and thus supporting the simpler form in Eq. (18), fitted above using the function 
/i(p^). Moreover, the x^/d-o.f. is not better for the 5-parameter fit compared to the 4-parameter fit, indicating that 
the latter is more stable. 

We thus conclude that our best fit is /i(p^) in Eq. (28), i.e. the 4d gluon-propagator lattice data favor the simplified 
expression in Eq. (18), implying (^^) = (ff)- 



VI. THE 3D CASE 

In this case the simplified fitting form fi{p^) in Eq. (28) is not able to describe well the lattice data. Indeed, even 
using improved momenta (see Table VII), the x^/d.o.f. values obtained are quite large. Moreover, as can be seen in 
Fig. 3, the fit clearly fails in the IR region.* The situation improves by considering the (5-parametcr) fitting function 
fiij)^) in Eq. (36) above, as can be seen from the results reported in Tables VIII and IX, obtained respectively using 
unimproved and improved momenta. Note that, as in the 4d case, the use of improved momenta helps to obtain a 
better fit to the data. 

One can also try to use the more general function 

U(v^) - C + + (42) 

- ^(p4 + ^2p2+i2)(p2 + ^)' (42) 

obtained by introducing the extra parameter I. In this case, the results of the fit using improved momenta for the 



Note that a fit using the more general form in Eq. (42) below is unstable in this case, yielding large errors for the fit parameters. 
Nevertheless, this fit suggests the factor (p^ + 1) in the numerator of (36), as adopted here. 

In order to highlight the results at small momenta, here and in Fig. 4 we present the plot with a logarithmic scale on both axes. 
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V 


C(GeV) 


M(GeV) 




s(GeV2) 


fc(GeV2) 


xVd.o.l 


140"' 


U.zoa (U.UUzj 


U.ooi (U.U22) 


0.552 (O.UUoj 


U.U18 (U.UU3) 


0.060 (U.UOd) 


10.48 


200^ 


0.289 (0.002) 


0.386 (0.019) 


0.552 (0.006) 


0.019 (0.003) 


0.032 (0.006) 


9.45 


240^ 


0.290 (0.002) 


0.393 (0.017) 


0.550 (0.005) 


0.020 (0.003) 


0.034 (0.005) 


6.55 


320* 


0.290 (0.002) 


0.389 (0.017) 


0.549 (0.005) 


0.019 (0.003) 


0.035 (0.005) 


2.89 



TABLE VIII. Fits of the gluon-propagator data in the 3d case, for different lattice volumes, using the fitting function /4(p^) 
in Eq. (36) and unimproved momenta [see Eq. (26)]. We report, besides the value of the fit parameters, the x^/d.o.f. obtained 
in each case. The whole range of momenta was considered for the fit. Errors shown in parentheses correspond to one standard 
deviation. 



V 


C(GeV) 


M(GeV) 






fc(GeV2) 


xVd.o.f. 


140* 


0.407 (0.001) 


0.654 (0.008) 


0.623 (0.004) 


0.022 (0.002) 


0.041 (0.003) 


2.14 


200* 


0.407 (0.001) 


0.655 (0.007) 


0.623 (0.004) 


0.024 (0.002) 


0.043 (0.003) 


1.92 


240* 


0.408 (0.001) 


0.662 (0.007) 


0.620 (0.004) 


0.025 (0.002) 


0.047 (0.003) 


1.59 


320* 


0.408 (0.001) 


0.656 (0.008) 


0.619 (0.005) 


0.023 (0.002) 


0.046 (0.004) 


1.19 



TABLE IX. Fits of the gluon-propagator data in the 3d case, for different lattice volumes, using the fitting function fnip^) in 
Eq. (36) and improved momenta [see Eq. (27)]. We report, besides the value of the fit parameters, the x^/d.o.f obtained in 
each case. The whole range of momenta was considered for the fit. Errors shown in parentheses correspond to one standard 
deviation. 



lattice volume V = 320^ are 

(43) 
(44) 
(45) 
(46) 
(47) 
(48) 

with x^/d.o.f. = 1.19. By noticing that I w 1 and by comparing these values to the corresponding ones from the fit 
using f4{p^) in Eq. (36) (reported in the last row of Table IX), it is clear that the two results are equivalent. At any 
rate, the values of x^/d.o.f. for the 6-paramctcr fit [using the function f^ij)"^) in Eq. (42)] and for the fit in Table IX 
are the same, indicating that the latter is more stable. 



c = 


0.405 ± 0.003 GcV 


u = 


0.692 ± 0.040 GeV 


t = 


0.635 ± 0.018 GeV^ 


s = 


0.025 ± 0.002 GeV^ 


k = 


0.050 ± 0.007 GeV^ 


1 = 


1.092 ± 0.103 GeV^ 



parameter 


propagation of error 


Monte Carlo analysis 


bootstrap analysis 


(GeV^) 


0.512 (1) 


0.512 (1) 


0.513 (1) 


(GcV^) 


-0.55(1) 


-0.55(1) 


-0.52(2) 


(GcV*) 


0.94 (1) 


0.94 (1) 


0.91 (3) 


pi (GeV") 


0.479 (2) 


0.479 (2) 


0.477 (2) 


P2 (GeV") 


0.09 (1) 


0.094 (9) 


0.100 (6) 



TABLE X. Estimates of the parameters of the RGZ propagator in Eq. (17) from fits (see last row of Table IX above) to the 

equivalent form /4(p") in Eq. (36). Errors arc obtained using propagation of error, a Monte Carlo analysis with 10000 samples 
and a bootstrap analysis with 500 samples. In all cases we considered the volume V = 320* and improved momenta. Errors 
shown in parentheses correspond to one standard deviation. 
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Q. 

Q 



0.1 



0.01 




FIG. 3. Plot of the 3d gluon propagator D{p^) (in GeV^^) as a function of the (improved) momentum p (in GeV) for the 
lattice volume V = 320^. We also show the fitting function /i(p^) [see Eq. (28)] with the parameters reported in the last row 
of Table VII. Note the logarithmic scale on both axes. 



In order to evaluate the condensates of the RGZ model, we thus consider only the results from the fit using /4(p^), 
given for the lattice size N = 320 in the last row of Table IX. (The corresponding plot is shown in Fig. 4.) By setting 
fiiv^) [see Eq. 36)] equal to the RGZ propagator (17) modulo the global factor C, we find (using propagation of 
error) the values for the condensates reported in the first column of Table X. Note that, using this fitting form, we 
arc able to evaluate M^, pi and \p\ (and therefore p2) separately. In this case, we can see that p2 7^ and p is indeed 
a complex quantity. This is consistent with the fact that the (four-parameter) fit to the simplified form fi{p^) fails, 
as seen above. 



Finally, we decompose the propagator as in Eq. (19) with (3 = a + ib, •y = a — ib, = v + iw e uj^ 



tw. I.e. 



parameter 


propagation of error 


Monte Carlo analysis 


bootstrap analysis 


a (GeV) 


-0.024(5) 


-0.024(5) 


-0.029(4) 


(Gc\--) 


0.0 f() (f) 


0.0 f() (f) 


().0i() (f) 


a (GeV) 


0.216 (3) 


0.216 (2) 


0.220 (4) 


&(GeV) 


0.27 (5) 


0.271 (3) 


0.275 (3) 


V (GeV^) 


0.215 (5) 


0.215 (5) 


0.23 (1) 


w (GeV^) 


0.580 (6) 


0.580 (6) 


0.57 (1) 



TABLE XI. Estimates of the parameters of the function fe{p^) [see Eq. (49)] from fits (see last row of Table IX above) to the 

equivalent form fiip^) in Eq. (36). Errors are obtained using propagation of error, a Monte Carlo analysis with 10000 samples 
and a bootstrap analysis with 500 samples. In all cases we considered the volume V = 320^ and improved momenta. Errors 
shown in parentheses correspond to one standard deviation. 
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Q. 

Q 



0.1 



0.01 




FIG. 4. Plot of the 3d gluon propagator D{p^) (in GeV^^) as a function of the (improved) momentum p (in GeV) for the 
lattice volume V = 320^. We also show the fitting function fiip^) [see Eq. (36)] with the parameters reported in the last row 
of Table IX. Note the logarithmic scale on both axes. 



we consider the function 



/6(P') 



2a + 2{av + bw) 



(49) 



We find (again using propagation of error) the results^ reported in the first column of Table XL Note that the 
imaginary part w of the complex-conjugate poles is more than twice the value of their real part v, as in the 4d case. 
Note also that the mass wi and the residue a associated with the real pole are very small. Moreover, a is negative, 
which may be associated with violation of reflection positivity, indicating that this mass cannot correspond to a 

physical degree of freedom. 

Also in this case we have repeated the analysis using a Monte Carlo estimate for the errors (with 10000 samples) 
and the bootstrap method (with 500 samples). The results are shown, respectively, in the second and in the third 
columns of Tables X and XI. The values of the fit parameters for the function (36) using the bootstrap method for 
the lattice volume V = 320^ are 



C = 0.411 ± 0.004 GeV 
u = 0.673 ±0.016 GeV 
t = 0.611 ± 0.006 GcV^ 
s = 0.025 ± 0.002 GeV^ 
k = 0.052 ± 0.005 GcV^ 



(50) 
(51) 
(52) 
(53) 
(54) 



which should be compared to the results shown in the last row of Table IX. Again, all results obtained using the three 
diff'erent analyses agree within one standard deviation. 



® Cleajrly, we have uil = k from fiip^). 
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We thus conclude that, in the 3d case, the data support a gluon propagator given by the general RGZ form (17), 
in which the condensate p is a complex quantity. This is in contrast with the 4d case seen in the previous section, for 
which p was real. There are also significant differences for the values of the other condensates and of in comparison 
with the 4d case. The masses from the complex-conjugate poles, on the contrary, have similar values in 3d and 4d 
(see Tables XI and IV respectively). 



VII. THE 2D CASE 

In the two-dimensional case the situation is different. Indeed, we know from Refs. [6, 56, 71] that the gluon 
propagator at zero momentum -D(O) does go to zero in the infinite-volume limit, even though one always has -D(O) > 
at finite lattice volume. Moreover, the behavior of D{p^) at small momenta is of the type p^, with some; non-integer 
power 1] « 0.8. This makes the fitting procedure more complicated than in the above cases, for which polynomial 
forms were used. After trying several generalizations of the fitting functions considered in the 4d and 3d cases, we 
found that a good fit to the gluon data can be obtained using the function 

which is a simple generalization of Eq. (28). Results of the fit for the various lattice volumes, using unimproved 
and improved momenta, are reported in Tables XII and XIII respectively.^*^ A plot of the fit for the lattice volume 

V = 320^ using improved momenta can be seen in Fig. 5. 

It is interesting to note that the function frip^) above can be decomposed as 



hyp ) = — 2 + — 2" = — ~ — ~ — 5 — ' (56) 



with 

a±{p'^) = a ± i{b + cp'^) , uj%=v±iw. (57) 

An estimate for these five parameters is reported in Table XIV, using again three different analyses for the error. The 
average values of the fit parameters in Eq. (55) using the bootstrap method for the lattice volume V = 320^ are 

C = 0.112 ± 0.001 GeV^ (58) 

u = 0.550 ± 0.013 GeV (59) 

t = 0.255 ± 0.006 GeV^ (60) 

s = 0.0152 ± 0.0008 GeV^ (61) 

Z = 0.326 ± 0.033 GeV^-" (62) 

7/ = 0.859 ± 0.026 . (63) 

We also tried an extrapolation to the infinite- volume limit of the gluon propagator at zero momentum D{0) using 
the function A + B/N" , where N is the lattice side in lattice units. This gives^^ 

^ = -0.002 ±0.010 (64) 

5 = 1.7 ±0.8 (65) 

zv = 0.7±0.1, (66) 



Note that in this case the use of improved momenta does not affect the quality of the fit significantly. Nevertheless, we choose to consider 
improved momenta in our analysis also in 2d. 

Recall that, with our notation, the 2d gluon propagator is dimensionless. 
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V 


C(GeV2) 


M(GeV) 


t(GeV2) 


s(GeV2) 


KGeV^-") 


n 


xVd.o.f. 


80^ 


0.073 (0.005) 


0.363 (0.041) 


0.265 (0.011) 


0.078 (0.008) 


0.403 (0.142) 


1.129 (0.151) 


2.92 


120^ 


0.069 (0.004) 


0.432 (0.027) 


0.252 (0.008) 


0.052 (0.005) 


0.566 (0.128) 


1.145 (0.088) 


2.75 


160^ 


0.067 (0.004) 


0.458 (0.022) 


0.250 (0.007) 


0.044 (0.003) 


0.665 (0.117) 


1.138 (0.066) 


2.68 


200^ 


0.068 (0.003) 


0.470 (0.022) 


0.254 (0.007) 


0.037 (0.003) 


O.o53 (0.104) 


1.088 (0.057) 


3.22 


240^ 


0.069 (0.002) 


0.469 (0.018) 


0.252 (0.005) 


0.031 (0.002) 


0.626 (0.076) 


1.051 (0.041) 


2.61 


280^ 


0.069 (0.002) 


0.483 (0.016) 


0.261 (0.005) 


0.029 (0.002) 


0.648 (0.064) 


0.994 (0.033) 


2.34 


320^ 


0.070 (0.002) 


0.483 (0.016) 


0.260 (0.005) 


0.025 (0.002) 


0.630 (0.062) 


0.981 (0.032) 


2.77 



TABLE XII. Fits of the gluon-propagator data in the 2d case, for different lattice volumes, using the fitting function frip^) in 
Eq. (55) and unimproved momenta [see Eq. (26)]. We report, besides the value of the fit parameters, the x'^/d.o.f. obtained in 
each case. The whole range of momenta was considered for the fit. Errors shown in parentheses correspond to one standard 
deviation. 



V 


C(GcV") 


u(GcV) 


t{GcY'^) 


s(GcV2) 


/(GcV^-") 


V 


xVd.o.f. 


80^ 


0.114 (0.002) 


0.433 (0.031) 


0.207 (0.012) 


0.031 (0.004) 


0.026 (0.043) 


0.684 (0.594) 


2.63 


120^ 


0.112 (0.003) 


0.486 (0.024) 


0.197 (0.008) 


0.020 (0.002) 


0.091 (0.050) 


1.003 (0.186) 


2.51 


160^ 


0.110 (0.002) 


0.503 (0.020) 


0.199 (0.006) 


0.018 (0.001) 


0.133 (0.044) 


1.027 (0.111) 


2.31 


200^ 


0.110 (0.002) 


0.523 (0.020) 


0.201 (0.007) 


0.015 (0.001) 


0.157 (0.046) 


1.017 (0.092) 


2.92 


240^ 


0.110 (0.002) 


0.519 (0.018) 


0.201 (0.006) 


0.013 (0.001) 


0.152 (0.038) 


0.955 (0.073) 


3.06 


280^ 


0.110 (0.002) 


0.530 (0.016) 


0.208 (0.006) 


0.012 (0.001) 


0.168 (0.033) 


0.904 (0.055) 


2.77 


320^ 


0.110 (0.001) 


0.539 (0.015) 


0.209 (0.006) 


0.011 (0.001) 


0.180 (0.033) 


0.909 (0.049) 


2.91 



TABLE XIII. Fits of the gluon-propagator data in the 2d case, for different lattice volumes, using the fitting function frip^) 

in Eq. (55) and improved momenta [see Eq. (27)]. We report, besides the value of the fit parameters, the x'^/d.o.f. obtained in 
each case. The whole range of momenta was considered for the fit. Errors shown in parentheses correspond to one standard 
deviation. 

with x^/d.o.f. = 1.07 (considering 7 data points). The fit improves if one sets A = 0. Indeed, in this case we find 

B = 1.9 ±0.2 (67) 
u = 0.71 ± 0.02 (68) 

with x^/d.o.f. = 0.87. Thus, the value of the parameter s in frip^) [see Eq. (55)] is consistent with zero at infinite 
volume, implying -D(O) = in the same limit. One should note, however, that the condition -D(O) = is not obtained 



parameter 


propagation of error 


Monte Carlo analysis 


bootstrap analysis 


a(GeV2) 


0.0550 (5) 


0.0550 (5) 


0.0559 (7) 


6(GeV2) 


-0.049(8) 


-0.049(7) 


-0.037(2) 




0.07 (1) 


0.07 (1) 


0.089 (8) 


V (GcV^) 


0.145 (8) 


0.145 (8) 


0.151 (7) 


w (GeV^) 


0.15 (2) 


0.15 (1) 


0.205 (6) 



TABLE XIV. Estimates of the parameters of the function fsip^) [see Eq. (56)] from fits (see last row of Table XIII above) 
to the equivalent form frip^) in Eq. (55). Errors are obtained using propagation of error, a Monte Carlo analysis with 10000 
samples and a bootstrap analysis with 500 samples. In all cases wc considered the volume V — 320^ and improved momenta. 
Errors shown in parentheses correspond to one standard deviation. Note that the value of rj can be obtained from the last row 
of Table XIII for propagation of error and Monte Carlo analysis, and from Eq. 63, for the bootstrap analysis. 
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FIG. 5. Plot of the 2d gluon propagator D{p^) as a function of the (improved) momentum p (in GeV) for the lattice volume 
V = 320^. We also show the fitting function frip^) [see Eq. (55)] with the parameters reported in the last row of Table XIII. 
Note that in 2d, with our convention, the gluon propagator D(p^) is dimensionless. Also note the logarithmic scale on the y 
axis. 



here with 6 = and ?; = [sec Eqs. (56) and (57)], i.e. with purely imaginary poles as in a Gribov-like propagator 
(12), but it seems to be due to the relations a = —b and v = w (see Table XIV). Thus, the behavior in the 2d case 
appears closer to the RGZ propagator (18) (with a non-integer power in the numerator) than to a GZ propagator (12). 
However, as mentioned in the Introduction, one cannot relate the fitting parameters to dimension-two condensates in 
the 2d case. 

Finally, taking into account the infinite-volume limit, we set s = in frip^)- In this case, the gluon propagator at 
small momenta behaves as D{p-^) w Clp^ /t^. Prom the conformal relation (see for example Ref. [56] with d=2) 

Dip") - (f,2)2K+(4-d)/2-l ^ (Q9^ 

we find the IR exponent k « 0.225 using jj w 0.9 (see Table XIII). This result is in reasonable agreement with the 
numerical estimate in Ref. [56]. 



VIII. CONCLUSIONS 

We have performed a systematic fitting analysis of Landau-gauge gluon-propagator data for SU(2) gauge theory in 

2, 3, and 4 space-time dimensions. The fit results were matched to analytic predictions from the RGZ framework, 
with the intent of calculating physical values for the dimension-two condensates in the theory, as well as to test the 
predictions in a neutral way. Indeed, as mentioned in the Introduction, our strategy has been firstly to find fits with 
convenient free parameters in simple (rational, except for 2d) forms and then to interpret the parameters in terms 
of physical quantities. In that way, we do not bias the fits by imposing relations between fit parameters from the 
predicted forms. In particular, a direct fit to the general RGZ formula [sec Eq. (17)] would involve six parameters 
(corresponding to the four condensates, the parameter A** and the overall normalization) with prescribed relations 
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among them. On the contrary, our best fits in the 4d and 3d cases involve, respectively, four and five free parameters 
(including the overall normalization), considering all generated data points. We find that the resulting fits in 3d 
and 4d agree remarkably well with the predictions from the RGZ scenario. 

In particular, the 4d results are well described by the simplified version of the RGZ gluon propagator in Eq. (18), 
equivalent to the simplest Gribov-Stingl form. This corresponds to a pair of complex-conjugate poles, as opposed to 
the Gribov propagator, in which the poles would be purely imaginary. Our fit results — using improved momenta 
and the fitting function /i(p^) in Eq. (28) — are given in Table II. The condensates and pole masses obtained for the 
largest lattice size are given respectively in Tables III and IV. The values for the condensates + pi , rm? and A* 
are in agreement with the ones obtained for the SU(3) case in Ref. [23], where the analysis was done using a slightly 
different notation, as previously explained. The quantitative agreement between the infrared limit of SU(2) and SU(3) 
theories was observed numerically before in [16, 75]. 

In 3d, our fits support the more general form of the RGZ propagator in Eq. (17). Fit results — using improved 
momenta and the fitting function f^ip^) in Eq. (36) — are given in Table IX, while condensates and pole masses from 
the largest lattice size are reported in Tables X and XI. In this case, the condensate p is a complex quantity and 
there are significant differences in the values of the other condensates and of A'' compared to the 4d case. Also, in 3d 
one has a real pole mass in addition to the pair of complex-conjugate poles. It is interesting to note that the masses 
from the complex poles assume similar values in 3d and 4d, with an imaginary part that is more than twice their real 
part. (We recall that a Gribov propagator would have a null real part.) Note also that the mass and the coefficient 
associated with the real pole in 3d are very small. 

In the 2d case, two particular features arise in the analysis of the data. Firstly, as known from previous lattice 
studies, the description of the infrared behavior of the gluon propagator D(p^) requires a non- integer power r/ of the 
momentum p. [See the fitting form f7{p^) in Eq. (55) and fit parameters obtained using improved momenta in Table 
XIII.] Secondly, the pole structure that best fits the data is similar to the one observed in the 3d and 4d cases — i.e. 
complex-conjugate poles with nonzero real part — at all considered (finite) lattice volumes. In the infinite-volume 
limit, one finds D(0) = 0, as would be the case for a Gribov propagator, with purely imaginary poles. However, in 
our case, the real part of the poles does not seem to vanish in this limit. The null value of -D(O) comes, instead, from 
an exact cancellation of the contributions from the two complex-conjugate poles (see values of pole masses in Table 
XIV). 

Our analysis strongly suggests — in d = 2, 3, 4 — a pole structure with complex-conjugate masses (with comparable 
rc;al and imaginary parts) for the infrared gluon propagator in Landau gauge. As stressed at the end of Section 
III A, one can interpret this result as describing an unstable particle [3, 35, 36]. In particular, by considering the 
position rrig — inigTg of the gluon pole, one can evaluate the gluon mass nig and its width Tg, which are in principle 
gauge-independent quantities [76, 77] . In our case, if we take as a reference the pole masses in the 4d case — i.e. the 
values uj% = v±iw with v w 0.3 GeV^ and w w 0.65 GeV^ from Table IV — we obtain 

mg w 550 MeV and Tg w 1180 MeV. (70) 

Note that the value for the gluon mass rUg is in agreement with other determinations [23, 28, 78, 79].^^ At the same 
time, the very large value for the width F^ would correspond to a lifetime Tg smaller than 10~^^ s, supporting the 
existence of very short-lived excitations of the gluon field. 

In summary, wc have presented fits — inspired by a renormalizablc action — allowing a good description of lattice 
data for the Landau-gauge SU(2) gluon propagator D{p'^). The data points range from about 4 GeV down to 20-40 
MeV, which are the smallest simulated momenta to date. Our results thus provide an accurate modeling of D{p'^) in 
the whole IR region, which will hopefully be a useful input in general studies of the IR sector of Yang-Mills theories. 



Wc recall, however, the importance of considering improved momenta [see Eq. (27)] for the fits. 
^•^ Let us mention that this complex-conjugate pole structure has been shown to describe also the longitudinal and transverse gluon 
propagators in Landau gauge at finite temperature [72-74] , at least up to twice the critical temperature . 

Let us remark that in Rcf. [80] the authors prefer to consider, instead of the gluon mass nig, the maximum wavelength of (confined) 
gluons, roughly corresponding to the inverse gluon mass 1 / rUg . 
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